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We propose a new method for determining total masses of low frequency eccentric binaries (such 
as, neutron star binaries with orbital frequency / > 10~^Hz) from their gravitational waves. In this 
method we use the frequency shift caused by periastron advance, and it works even at low frequency 
band where chirp signal due to radiation reaction is difficult to be measured. It is shown that the 
total masses of several Galactic neutron star binaries might be measured accurately (within a few 
percent error) by LISA with operation period of ~ 10 years. 



O ■ I- INTRODUCTION 

o 



The Galactic binaries such as close white dwarf binarie s (CWDBs) or neutron star binaries (NBs) are promising 



^ . targets of the Laser Interferometer Space Antenna (LISA, http:/ /hsa.jpl.nasa.go-y , [|1|). One of the important aim of 
) ^ ■ gravitational-wave astronomy is to extract out information of sources that emit gravitational radiation. In the case of 
, a Galactic binary we want to know masses of two stars, orbital parameters (semi- major axis, eccentricity, inclination), 
^\ • distance to the binary and so on. At the final in-spiral phase {e.g. the last three minutes of NBs) we can, in principle, 
OA . determine various parameters by fitting time evolution of wave signal with post-Newtonian expansion Q. But at 
LISA band many compact binaries evolve very slowly and situation is largely different. 
' It is observationally known that NBs (and some NS-WD binaries) with orbital period ^ Iday have large eccentricities 
J — 1^. This is explained by the "kick" effect and instantaneous mass loss at birth of a neutron star, in contrast to CWDBs 
whose orbits are circulaized by strong tidal interaction during mass transfer phases. As the eccentricity e decreases 
with increase of the orbital frequency due to radiation reaction Q, its effect is supposed to be negligible for most 
^ , high-frequency sources that will be searched by ground-based detectors (TAMA300, GEO600, LIGO and VIRGO). 
' But the eccentricity of NBs can be e ~ 0.1 at frequency / ^ lO^'^Hz as expected from PSR B1913+16 and it 
I might cause observable effects for LISA. 

' Frequency series of gravitational wave from an eccentric binary is affected by the periastron advance |^ (see Ref. 
'■q^' for radio observation of binary pulsars) . Using this fact we might obtain information of binaries even at low frequency 
I band where chirp signal would be difficult to be measured. In this Letter we propose a new method for determining 
O \ binary parameters (mainly total mass), and investigate its prospect for studying Galactic NBs with LISA. 
• 

■ 

II. GRAVITATIONAL WAVE FROM ELLIPTICAL ORBIT 

' Let us study gravitational waves from elliptical orbits following Ref. Q (see also Q ) . With quadrupole formula of 
^ I gravitational radiation ||^ two polarization modes in TT-gauge are written as follows 
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where O represents direction of the orbital angular momentum (inclination) and $ represents direction of the pe- 
riastron in the orbital plane (we put $ = when the periastron is on the plane determined by two vectors: the 
orbital angular momentum vector and the direction vector to the observer ||^). The amplitude /iq is given as 
/lo = 4(27r)~^/^G^/'^c~''m^{f^pr~^(Pb)~^/'^ {rrichirp. chirp mass, r: distance to the binary, Pb: orbital period from 
periastron to periastron). We have defined the frequency series /„ = Pj^^n. The angle changes due to the pe- 
riastron advance and we can denote it as 2$ = 5/ x < by adjusting origin of the time coordinate t. Thus the 
periastron advance causes frequency shift 5f that is given in terms of the total mass rritotai and the eccentricity e as 
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5f = 6(27r)^/^(Pf,)~^/^G'^/'^c~^mj^j^^j(l — e^)~^ (twice of the frequency for the periastron advance |^). The frequency 
fn now sphts into a triplet (/„ — 6f, fn, fn + Sf). Effect of the angle 2<I> is related to rotation of a coordinate system 
and appears in the same form for every n-mode. This is an essential point. The simple replacement 2$ 26 f x t 
into expressions (|l|) and (||) corresponds to an approximation that neglects terms of 0{PbSf) <C 1 |^. The coefficients 
Sm Cn and An are given by the n-th Bessel function Jn{x) and its derivative J„{x) = dxJn{x) as 

2(1-6^)1/^ , 2n(l-e2)^/2 
A„ = J„(ne), Sn = '^n('^e)H 5 Jn{ne), (3) 

Cn = — Jn {ne) + ^ 1 J„ (ne) . (4) 

e 

We can expand these coefficients around e = (circular orbit) and find that only n < 3 modes have terms 0(e°) 
or O(e^). They are given as Si = -|e + O(e^), Ci = -|e + 0(6^), Ai = + ©(e^) for n 1-mode, S'2 = 
1 - |e2 + 0(e4), C2 = 1 - fe^ + 0(6*), A2 = ^e^ + 0(e4) for n = 2-mode, and 5*3 = -|e + 0(e3), C3 = 
— |e + 0(6"^), ^3 = 0(e'') for n = 3-mode. Our basic strategy is to compare 0(e°)-term of n = 2-mode and 0{e^)- 
term of n = 3-mode (as we see below, the eccentricity relevant for our analysis is expected to be small e ^ 0.1). From 
these two frequencies /2 — 5f and fs, — Sf we obtain 

3(/2 - Sf)/2 - (/3 - 6f) - -5f/2, (5) 

and thus the total mass mtotai is estimated for a binary with small eccentricity e (note also we might determine the 
eccentricity e by comparing amplitudes of two waves). Amplitude of n = 1-mode is smaller than n — 3-mode, and 
furthermore the binary confusion noise would be larger at lower frequency Therefore we investigate prospect of 
our method by studying detectability of n = 3-mode in stead of n = 1-mode. 



III. DETECTABILITY 

The frequency difference 5f due to periastron advance is expressed as j|] 

1.2x10-^ /m,,,,A^/7 h \"\ 
^^^^^[jm-J UxlO-^Hz j 

and the accumulated frequency difference (chirp signal) due to gravitational radiation reaction within observational 
period Tobs is given as [Q 

■ 3.0 X 10-« ( m,mrp \'''^ ( /2 ^ A , 73 2 , 37 ,\ ( T,bs\ ^ 
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Thus at lower frequency the difference due to the periastron advance can take larger value. The estimation error 
(resolution) A/ for wave frequency / in matched filtering analysis is given by the Fisher information matrix of fitting 
parameters. For the set of three unknown parameters (/, /, initial phase) of quasi-monochromatic wave (/Tohs <C /) 
we obtain (§ [] 

A/ = AVi.-'T2SNR-' = 6.6 X IQ-^" Hz, (8) 

where SNR is signal to noise ratio of the wave. For the secular frequency modulation / we have estimation error 
A/ = ?,\/b'K~^T~^^SN R^^ . From equations (|^)(^) and (||) the frequency shift 5f would be resolved 
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'^When we include the direction and orientation of the source in fitting parameters, the estimation errors might be somewhat 
larger than our evaluation. We should also notice that these angular parameters are fixed well by / = 2-mode and we might 
not need to fit them for determining frequency of / = 3 mode. 
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Here (and hereafter) we consider contribution of the estimation error for 5f only from n = 3-mode, as SNR of 
n = 2-mode would be larger. In the foUowings we estimate the number of Galactic NBs whose n ~ 3-mode is detected 
with SNR > 10. 

First we discuss time evolution of the eccentricity parameter e. Using quadrupole formula for gravita- 
tional radiation the orbital semi-major radius a is related to the eccentricity as a/ui = (1 — ef)/(l — 

e^) (e/e,)^^^^^ [(1 + 121eV304)/(l 121e2/304)]**™/^^^^ where a, and a are their initial values [|. At a/ai < 1 

we have an useful approximation e ^ (a/ai(l — ))^^^^^ = [f / fi) ^^^^^ (1 — ef)^^^/^^ei where Kepler law is used. 
For a given orbital period the eccentricity should have some distribution function F{e\ Pb) [|o|, but it is poorly known 
at present, as only several NBs have detected so far Here we use the binary pulsar PSR B1913-I-16 as a refer- 
ence. This system has orbital period Pb = 2.8 x lO^'sec and eccentricity e = 0.62 ||]. From above approximation its 
eccentricity evolves as 

e ~ 0.13(/3/0.001Hz)-i9/i^ (10) 

Hereafter we use this 6-/3 relation in our order estimation. 

Next we study the frequency and spatial distribution of the Galactic NBs. Assuming that the former is in steady 
state, we can evaluate the distribution function dN/df using the coalescence rate Rns of Galactic NBs as 

dN ^ fdfY' 00, ,,^4^ Rns \ f f2 A -1 



Kalogera et al. ||ll[] recently estimated the event rate as Rns — 10^® — 5 x lO^^yr^^ (see also [Q). For spatial 
distribution of Galactic NBs we use the standard exponential disk model p{R, z) — p^exp {—R/ Rq) exp {—\z\/ zq) , 
where {R,z) is the Galactic cylindrical coordinate |]ll|,^. We fix the radial scale length Rq = 3.5kpc and the disk 
scale height zq = 500pc, and assume that the solar system exists at the position R = 8.5kpc and \z\ = 30pc. 

We have a relation h^ifz) — — 9e/i2(2/3/3)/4 between n — 2-mode and n — 3-mode (see eqs.[||[|]). The amplitude 
oi n — 3-mode is given as 

where we have taken angular average with respect to orientation of sources . From equations ([l^) and ( [l^ ) we can 
estimate the distance to a binary whose n — 3-mode is detected with a given SNR, 

lonnf f^ \ ''^^^ fhrms{f?,-,Tobs)\ ^ ( \ n o^ 

^^^'"% 3xlO-3Hz j y 10-^3 j (^^jP^' (13) 

where hrms{f3,Tobs) is the noise spectrum. We use the angular averaged sensitivity (effectively a factor of a/5 
degradation 10,0) to take rotation of LISA into account, and adopt the noise spectrum given in Ref. |l5|. When 
the effective frequency bin TJ^J is occupied by more than one Galactic CWDBs (more precisely, number of fitting 
parameters is larger than that of data), their confusion noise becomes important and the sensitivity hrmsif,Tobs) 
would be significantly worse at lower frequency. Then binaries very close to the solar system would be only resolved. 
This transition occurs at frequency ft — 1.6 x 10"'^ (T/10yr)~'^^^^Hz for abundance of Galactic CWDBs estimated 
in Ref. [l5|] . As we see below, the prospect of our method depends strongly on the transition frequency and the 
position of ft is very important. In an extreme model that the Galactic halo MACHOs are WDs, the frequency ft 
can be close to ~ 10~^Hz |16 (see also fl^). If this is true, our method would be severely limited. Spatial filters that 



depends on the angular position of sources might work effectively in some cases |18|| . At f ^ ft the noise spectrum 
hrmsif, Tabs) is mainly determined by the detector's intrinsic noise and the confusion noise by extra-Galactic CWDBs 
that depends on cosmological evolution of binary systems |H,Ojl9] . At this frequency region the total noise spectrum 



behaves simply as h„ns(.f ,Tobs) oc T^^^J . In the present analysis we only count NBs with frequency /2 = 2/3/3 > ft- 
Now we can estimate the number of NBs whose n = 3-mode is detectable by LISA with SNR > 10. In figure 1 
we show the results for observational period Tobs — lyr and Tabs = lOyr. The total number (integrated in frequency 
space) becomes 0.11(i?Ars/10~^yr) for lyr observation and 3.7(i?Ars/10~^yr) for lOyr. Thus our method would be 
effective for a long (but realistic) observational period Q. For these NBs {Tobs = lOyr, SNR < 10 and /2 > ft) 
the total masses nitotai would be resolved within a few percent accuracy (see eq.|^). In figure 1 we also plot all the 
Galactic NBs (dashed-line). Note that our method works well for Galactic NBs with f > ft in the case of Tobs = lOyr. 
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The estimation error for the overall amplitude of wave becomes ^ 0.15{SNR/10)~^ in the LISA band (see table 4.5 
of [0). Thus the eccentricity parameter e for a detected NB (with SNR > 10 for n = 3-mode) would be measured 
within ^ 15% accuracy. It has been observationally clarified that the mass of a neutron star is ~ 1.4M0 (that 
of an ordinal white dwarf is ^ 1.2Mq). If most of an eccentric binaries with mtotai ^ 2.8Mq and / ^ 10~'^IIz are 
either NS+NS or NS-WD binary we can obtain further information of such binaries detected by our method in 
the following manner. One of the two stars would be a neutron star with mass ~ 1.4Mq, and we can estimate the 
chirp mass of the system. Then the distance r to the binary is obtained as proposed by Schutz |^ . 

Using a similar method based on the expression A/ given just after equation (||), we have also estimated the 
number of Galactic NBs whose chirp signal few for n = 2-mode due to gravitational radiation are measured within 
5% accuracy (5 x 3/5 ~ 3% for the chirp mass eq. 01) [||. As shown in figure 1, all Galactic NBs with /2 > 4 x lO'^Hz 
satisfies this observational criteria for Tobs = lOyr. There might be ~ l(i?Ars/10~^yr~^) NBs whose individual masses 
can be determined accurately from observed mtotai and rrichirp- 

One might wonder whether two gravitational waves from different binaries (mainly CWDBs) are confused as n = 2- 
and n = 3-modes of a same eccentric source. Frequency distribution of Galactic close white dwarf binaries (CWDBs) 
is estimated as dN/df ^ 1.4 x 10^(//10~'^IIz)~^^/'^IIz"^ |15|. We can specify the direction of a monochromatic source 
with estimation error < 10~^sr for frequency / > 2x lO^'^Hz (in the case of signal to noise ratio:10, table 4.5 of Ref. 0|) 
using annual modulation of gravitational wave due to motion of LISA |l|,|2^. The mean frequency interval {Sf)int for 
binaries within a box of - lO^^sr is given as (5/)™* - {IQ-^ /An)'^ {dN/df )-'^ - 6 x 10-^{f/A x 10-3Hz)"/3Hz. The 
typical frequency difference due to periastron advance df (eq.^]) is much smaller than this interval and coincidence 
of miss-identification would be {Sf)/{Sf)int ~ 2 x IQ-^ifa/S x lO-^Hz)-^ The total number of Galactic CWDBs 
within frequency 2 x lO^'^Hz < /2 < 4 x 10^'^Hz is estimated as ^ 4 x 10*. Thus number of miss-identified CWDB 
pairs would be at most ~ 80. The estimation error of LISA for the direction of the orbital angular momentum is 
- IQ-isr {SNR = 10, table 4.5 of Ref. |]). Thus there would be only ~ 80 x 10-V47r - 1 miss-identified CWDBs. 
It would be also possible to distinguish NBs using the estimated total mass itself, as most CWDBs have total masses 
smaller than 2.4M0 ||l|l. 

Let us summarize this Letter. We have proposed a new method to determine the total masses of low frequency 
eccentric binaries form their gravitational waves. We use effects of periastron advance on the frequency space of 
gravitational wave. Our method would be effective for Galactic NBs (/ > 2 x 10~^IIz) with a long term (^ lOyr) 
operation of LISA. 
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FIG. 1. Various distributions of Galactic NBs within frequency bins /o < /2 < 1.26/o. The dashed-line represents all 
Galactic NBs. The filled triangles are the number of NBs whose gravitational wave of n = 3-mode are detected with SNR > 10, 
and open triangles are those whose chirp signals {f)Gw due to radiation reaction are measured within 5% accuracy. We take 



the coalescence rate of Galactic NBs at 
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